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I. INTRODUCTION 


A. THE PROBLEM 

The problem considered here is the search for a randomly moving target with a 
fixed speed and direction. We assume that the initial target distribution is uniform 
over a very large search area, and the searcher, carrving a cookie cutter sensor, is 
conducting a systematic search Within that area. The question addressed here is, what 
is the expected area swept, A(t), after t hours of search. We need to express our 
answer in terms of known parameters of target speed (u), searcher speed (v), sensor 
radius (R), and search time (t). 

A standard approach to this problem is to assume that the expected relative 


speed between a searcher at speed v and a target at speed u 1s 


5 


0 TE Re 
w = (1,27) J, v (u- + ¥-- 2uvcos@) dé. 
Tien the expected area sedreied by time) tis A(t) = 2Rwt, where Roismaie 


detection range of the searcher. This model was introduced by Koopman in Search 
and Screening [Ref. 1]. In this thesis we will integrate (1-p) over the area to calculate 
the equivalent area searched (A,), where p is the posterior target density for any point 
in the area. The method 1s to derive analytic expressions for the equivalent area swept 
(when possible), to generate target density curves for the area, and then simulate some 
of the main results. 
During this study, the following assumptions hold for each target at anv time t. 
[. Initial target position 1s uniforniy distributed over the ocean. 


2. Target course is uniformly disributed over (Oi): 


(2 


. Target position and course are independent random variables. 
Wy 


stationary searcher and stationary target cases. In addition, when both target and 


Om 


model used in this project gives results agreeing with Koopman’s model in 


searcher are moving with fixed speeds, the model developed here shows that the 
equivalent area swept is, contrary to intuition, the same for all v S u in a specified 
area centered about the searcher. Moreover, the specific posterior target density 
profile is computed and may be of use in planning operations such as follow-on 


searches. 


Il. CALCULATION OF TOTAL DEPLETED AREA 


A. THE MODEL 

When a searcher is progressing on its course at the constant velocity v among a 
uniform random distribution of targets of speed u, it is frequently important to know 
the probability of a particular target passing within the stated range of R mules of the 
observer and hence being detected. In some cases R may be the range within which 
the searcher can sight the target. [n others, it might be the range within which the 
target can detect the searcher’s presence. Again R may be effective gun fire range of 
ie searcher against the target, or vice versa. he problem we consider is this : Given 
a target located at (r.B), what is the probability P(r.) that the target in the past came 
within range R of the searcher? 

Evidently P depends on (r.B) and P=1 if rSR. When r2R, the target has 
entered the circle if and onlv if its relative velocity vector Ww points to the detection 
circle. The angular range of the target course giving a detection sometime in the past 
is drawn in Fig. 2.1 (shaded angle). Then, given the assumption that target course 1s 
uniformly distributed between 0 and 2m radians, P(r.B) becomes the ratio of the angular 
fame@ceo! U to: 2m. The problem is thus reduced to the geometry of Fig. 2.1. and the 
formuta for P(r.B) is obtained by straightforward trigonometry. 

The expression of P(r.B) is as fottows: P(r.B)= O(r.B). 2m where O(r.B) is the total 
range in radians of u. This method was used by Koopman in Search and Screening 
[Ref. 1] Koopman calculated the probability P(r.8) that the target at (r, §) will in the 
future come within range R of the searcher. The geometry in Koopman’s problem is 
very sinular to that studied here. In both cases a detection occurs whenever an 
extension of the target relative velocity vector comes within range R of the searcher. 
In Koopmian’s case, the relative velocity vector is extended forward to intersect the 
detection disk. Here it is extended backwards. 

For tater comparison with our more general result we will first took at the two 
Scinemie Cases in which either the target or the searcher is stationarv. We will be 


assuming a prior target density, p. of | target per unit area. 


; RW w 
VY 


Figure 2.1 The Probability of Entering Circle. 





B. STATIONARY TARGET 

As seen in Fig. 2.2, if the target is fixed, the searcher develops a rectangular 
“target-free tail” of area 2Rvt, where t is search time. In this case target density p is 0 
inside the target free tail and on the disk, and it is 1 outside these areas. The total 


swept area at time t, A,{t), 1s mR* + 2Rvt. 


C. STATIONARY SEARCHER 
Now let us consider the case in which the searcher 1s fixed. For a point at the 
distance r S ut+R from the center of the detection disk, the posterior target density 1s 
reduced by an amount proportional to the shaded angle in Fig. 2.1 Hence, 
U8 ame — 


Now simple geometry yields, 


ArCSimuns 1 
; es ) , rs ¥ (ut)?+R? 
1 


p(r) = 
l Nes Bien 





Figure 2.2 Searcher Motion In Stationary Target Area. 


Note that we have not given an expression for p(r) when r 1s in the interval 
(/ (ut)?+R* , ut+R). Here the detection disk “edge effects” determine p(r), as 
suggested in Fig. 2.3. For mathematical convenience we will generally ignore these 
edge effects and use the arcsine formula above to calculate p(r) for all r<(utt+r). The 
seeeeciicct will be small for ut>=>R. Note that p(r) is 1 for all r>utt+R. This is 
because no targets can be captured by time t that start at a greater distance than 


ieut. INeglecting the edge effect we have, 





loune 25 


A(t) = Fe Cpl) rdrd@ + mR? 


on R+ 
ay): = (in)f™ i Utarcsin(R/r) rdrd@ + mR- 


_ 
ae Ie Marcsin(R/r) rdr +mR? (eqn 2p 


After integrating the Equ. 2.1 (see [Ref. 2] ) we substitute arcsin( Ri) 


(good when r>>R) and get 
A(t) = uRt{l+ /(1+2Ryut)} + R21 + 2/2) 


Eqn 2.1 overestimates the actual swept area because of the above mentioned edge 
effects on the detection disk. A lower bound can be found by evaluating the integral in 
Eqn. 2.1 only between R and ./ R*+(ut)*. Fig. 2.4 shows a plot of these upper and 
lower bounds, together with a plot of 2Rut+mR* which, according to Koopman’s 


result mentioned in Chapter I, should be the actual searched area by time t. 


§ 
° 
i) 
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Figure 2.4 


It is interesting to take the derivative of Eqn. 2.1 with respect to t to calculate 


the sweep rate. 


dA (t) “ 


A 20 (Reet) arcsin( R/(R + ut)), 
t 


which for small R and large t becomes 





dA (t) R 
——— _ ~ 2u(uttr) = 2uR 
dt (ele ee) 
Fience, 
A(t) =2uRt (eqn 


We note that using this methodology when either the searcher or target is 
stationary is not particularly recommended, but is presented here to show that the 
methodology gives reasonable results in these simple cases and to demonstrate the 


method for a case in which the correct answer 1s known. 


D. TARGET AND SEARCHER BOTH MOVING 

In this problem both the target and the searcher are moving (Fig. 2.5). @ is the 
angle between the horizontal line and the target velocity vector u (corresponding to 8, 
or @, in Fig. 2.1). 





Figure 2.5 Relative Velocity Vector Diagram. 


In the frame moving along with the searcher, targets starting from any given 
point can move in any arbitrary direction, but only the ones with velocity vector u 


inside of the shaded area in Fig. 2.1 could have been captured. So we need to find 6. 


il 


in the moving frame, and transform it to 9. in the stationary frame, to find the angle 


shielded bv the disk for each point im the area. in Tircameres 
-W> V=wvcos@ (eqn 2.3) 


If W=U-V is substituted in Equ. 2.3 we get 
-(U-V) Vo = wy cos@’ 
VeV-UV = wv cos (eqn 2.4) 


a 


Since V° V = vw~ and U* V = -uvcos@, substituting into Equ. 2.4 yields 


; 
v~- + uvcos8 = wvcos8 


. (v+u cos®) : 
ee ————— (eqn 2.5) 
Ww 


By using the Cosine Law we get 


ty 


2 ae 
Wet Us =e case 


5 > 
= ur + vr +2uvcos8 


a . : . = ° ° . y 
Now squaring both sides of the Eqn. 2.5 and substituting for w-, 


> > + 
(¥> 4 U- Cos Ge eunicesu) 


cos" = (eqn 2.6) 


a | =. 
(Os > scoce! 
? ayn 5 7. 3 @ 2 
u“cos-8 + v-cos*@ + 2uvcosGecoss@ -v- = Uccos-G - 2Uuvcosus— a" 


5 9, nN e > + ay. 2 
u-~cos-§ + 2uvsin-8 cos - (u-+v~)cos-@ + v~ = 0 


The solution of this quadratic equation gives, 





OOS ae f-uvsin7@ + /(uvsin7@ )* + u*(cos70 (u” + v7)-v7)} / ur (eqn 2.7) 


-_ 


ee 


So when we have @ in the moving frame, we can transform it back to 9, in 
stationary frame. @ can be found by trigonometry. For the following calculation let 
us take L as the horizontal axis, X 1s the vertical axis and the bottom of the detection 


circle is the coordinate center (Fig. 2.6). Using the notation in Fig 2.6, 


Figure 2.6 


Sina@=R// and tanB = R-x/ Y. 

Since /? = (R-x)? + Y?, 
@=arcsin{r / J (R-x)?+ ¥2} and B= arctan {(R-x)/Y} 
0 =a+B and -0, =a-B 


By subtituting the @ and B values in this formula we get 


0,’ = a+ B = arcsin(R / /(R-x)°+Y7} + arctan{| R-x| / Y} 
-6,° = a-B = arcsin{R / /(R-x)*+ Y7} - arctan{|R-x| / Y} (eqn 2.8) 


rd 


So we have calculated @ 1 >: We can now transform these angles to 0, , by using 
our transformation formula, Eqn. 2.7. The shielded angle for anv given point in the 


area will be |8, - 8,|, and target density for that point is p = {1- |@,-0,| / 27}. 


13 


Earlier, we selected a point on one side of the geometric tail and derived Equ. 2.8. 
For symmetric points on the opposite side of the tail, 0, and 0’, do not change (Fig. 
2.7) and Baie 2-o) sulllemolds: 





Figure 237 





Figure 2s 


If we consider the same problem for the region outside of the geometric tail (Fig. 2.8) 


and derive an analvtic expression for 0’, and 0", we hae 


= arcsin{[R / /(R-x)’ + Y*] and 6B = arctan [(R-x) / Y] 
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6, = a + B = arcsin{[R . J(R-x) + Y7] + arctan{(R-x) Y] 


{= @-B = arcsinf[R /(R-x)* + Y"] - arctan{(R-x) | Y] 


These are the same expressions obtained in Equ. 2.8. So these equations can be used 


aes) 
l 


to find the density p = (1-|0, - 8,| / 2m} for any given point. 

We cannot find the total effective area searched (A.) by closed form integration 
as we did before when the searcher was fixed. Instead we find 9’ and transform it to 8, 
then calculate total effective area searched by numerical integration. Let k=v,u for 
computational purposes. 

We can use the computer programis given in Appendix B to calculate densities for 
each point in the area considered. Plotting of these values gives density curves in Fig. 
2,9 for the case in which target is moving twice as fast as the searcher. Fig. 2.9 and 
2.11 show that density curves propagate lke a sound or shock wave in the air. Note 
that the horizontal and vertical scales are different in these figures. The contours in 
ice. 7 aresocw | walmcicseminmi ie 2th i we look closely at the forward-most 
contours and the target free tail (the area in which p=0), we can see that they follow 
the Cherenkov formula [Ref. 3] for shock front directions. We will see the shape of 
the target free tail in more detail in Appendix A. 

When k= 2, the searcher is moving twice as fast as the target. and the searcher 
develops an empty tail behind. This target free tail starts to develop for k> 1 and gets 
larger and larger as k increases as 1s seen in Fig. 2.13. 

It is important to note that the target densities calculated are steady state values. 
That is, they would result only after the search has been in progress for a sufficiently 
long time. A “sufficiently long time” for any point x 1s the maximum time required for 
a target starting on the edge of the detection disk to reach point x in searcher- 
stauonary relative space. 

After running the programs in Appendix B for u=2 and v=1 and plotting them, 
we obtain the graphs given in Fig. 2.9 and 2.10. Now the question is what is the total 
swept area A? The total swept area after search for a time t (where ut> > R. so we 
can neglect edge effects) is obtained by integrating (l-p) over the entire search area. 
Since this is not possible to do analytically, we use a numerical integration method. 
According to Koopman’s model, this area will be proportional to w. Our results and 


the simulation presented in Appendix C are consistent with this. 
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Figure 2.10 Density Surface k=.5. 
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Figure 2.13 Development of the Tail. 
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However, we have found that integrating 1-p over any finite area centered about the 
presemt searcher position yields aeconstant for any vS u, for a given u. This results 
because the amount of increase in density in front of the searcher is the same as the 
PmicunurOrecccrcase InTGeNsitymbcmind the searcher for a fixed target speed u 2 V. 
That results in the same total effective area searched whenever the relevant area A has 
a front-to-back svmmetry about the searcher. 

To see why this can be so, in Figure 2.14 select one point in front of the searcher 
(K°) and another point behind the searcher (K) which 1s symmetric with respect to the 
searcher. Now take a close look at these points to see what is happening as v changes 


and u remains fixed. 





Figure 2.14 


In Fig. 2.14a we draw a velocity vector diagram for v=1 and u=3 to see @ and 
Q’ angles which we had calculated by computer program. @ and F represent the 
Epicldea sector by the searcher. Now let us keep u=3 and change v from I to 2. As 
seen in Fig. 2.14b, as v increases, 8 increases but F decreases. In order for (F + 98) to 
remain constant, the amount of increase behind the searcher must be equal to the 
amount of decrease in front of the searcher. Hence, the increase in p at the point 
behind the searcher is compensated by a decrease in p at the point ahead of the 
searcher. This gives us an intuitive idea about what is occuring as v increases and u 1s 


Kept fixed, but we still need to prove it. 


| 


In Section II.d we picked a point K (Fig. 2.14) behind the searcher and derived a 


formula for cosO (eqn 2.7). By using that formula and letting k= v u we get 
pe 
0) 5 = arccos{ k{ -sin-8’ + cos®’ // (k*-sin-6’) ] } {eqn 2.9) 
Now let us apply the same method to the point K’ in front of the searcher and 
derive an expression for F. In Fig. 2.14, 
W = U-V 
-W: V = wvcos(tt - 8’), so 


W: V = wvcos@ 


Also, 
We = (U-V v= aay (eqn Zoi 


Substitute VV = v* And U: V = uvcosF in CGT tO 10 eld 


ee ) 
We\ = euxcose —suvecosic=..- 
ucosF-v 
cos0’ = (eqn 2.11) 
Ww 


: > . @ 
By Using cOSine theorem. WwW) = Ur cee =u cos 


2, a 2 
ey u*cos”F + v4#-2uvcosk 
cos-Q° = a 
Vs ty cet Coss 


Solving this equation for cosF gives 
cosF, » = k[sin°®’ + cos®/(k>-sin-@’)] 
Eee tecos ak! sin-8’ + cos //( k--sin-8’) ] }. 


If we look at the argument inside of the { } in this expression, we notice that it 1s 
the negative of the argument inside of the { }simgequiezass 


Now applying the chain law, 


tJ 
tw 











d6 d d . SESE 
a = —— arccos(x) ty,u[ -sin?@’ + cosO’ / (v-u)*-sin-@' | } 
dv dx ug 
d Ure ; > 5. MAN 
= - —— arccos(x) fy u[ sin-@’ = cosO’ // (v.u)--sin-8 } 
dx Y 
7 CU tig 
ale 


That means, as searcher speed tncreases, the amount of tncrease of shrelded angle 
behind the target ts equal to the amount of decrease of shtelded angle in front of the 
searcher. Since the amount of change of shielded angle at each pomt will be 
compensated by s¥mimetrical points with respect to the searcher and swept area 1s 
equal to the sum of shtelded angles over the area, there will be no change tn equivalent 
area swept for a symmetrical area around the searcher’s present position. This result 
does not mean that we will get the same total equivalent swept area for different vSu 
values at a particular trme t. The calculation of that area tmtegrates (1-p) over a 
circular area of radrus R+ut centered on the initial searcher position. If we consider 
the time dimension and the searchers inittal positron, the equivalent swept area will in 


general change with u, Vv and t. 


ty 
oe) 


iil. CONCEUSIGN 


In searching in an area, a major concern is to calculate the swept area by the 
time t. This information can be used to determine how effective the search is, the 
number of searchers to be allocated to this area and how tong this search should be 
conducted. Koopman approached this problem with the expression for swept area. 
A(t) = 2Rwt, where w is the expected relative speed. In this project we exaniineauume 
same model, with a different approach calculating the total swept area as a function of 
target and searcher parameters. The reason we used this model ts that the area swept 
Is not just that behind the searcher unless the target is stationary. We need to know 
how this search changes the initial target densitv at different points inside of the 
affected area. 

Since we considered a very large ocean, we did not attempt to computememe 
percentage of the area swept. We defined the MOE as either the total area swept by 
the time t or equivalently, the number of detections during that time. assuming a 
uniform initial target densitv. For extreme cases (stationary target and stationary 
searcher) we found the total swept area analytically as a function of the relevant 
paramieters obtaining the same total swept area as in the Koopman model (which 
calculates the total swept area bv using expected relative speed). Suprisingly, we saw 
that when the target and the searcher are both moving, the amount of change in 
densitv in front of the searcher is the same as the amount of change in density behind 
the searcher, as v changes for a fixed u2v. For a search area, with front-to-o7e. 
svmiumetry about the searcher’s present psition, and after a sufficiently long search 
(defined in Chapter II ), the total swept area is constant for any search speed v less 
than or equal to the target speed u. 

The analvtic result for total swept area bv the time t when doth the target and 
the searcher Were moving was confirmed by use of the simulation program in 
Appendix C and is consistent with our detailed numerical evaluation of posterior target 
densities. Also the target density curves can be useful for deternuning optimal search 
patterns and efficient locations for the other searchers. Other extensions of our 
analvsis include targets moving with an arbitrary distribution of speed. Through 


further analysis using an evaluation of posterior target densities this distribution of 


target speed could determine the optimal search speed and search pattern. These 


applications have not been pursued here. 


APPENDIX A 
THE SHAPE AND THE AREA OF THE TARGET FREE TAIL 


Let us examine the target free tail that the searcher develops behind. If we look 
at our transformation formula (eqn 2.7), we can see that this transformation will not 
work inside the target free tail. This is because any point inside this tail is shielded in 
all directions. In other words there is no way for any of the targets in the area to get 
in this tail. This happens only when the discriminant of equation 2.7 is negative. So if 
we make this discriminant equal to zero, we can determine the boundaries of the target 


free tail. 
(uvsin?@’)? + u*(cos*@(U2+ v2) - v2) = 0 
v> sint@’ - (u? + y*)sin70 Son) 
By solving this equation we get 
sind, = u/v and sin®,’ = -u/v 
This the same formula used to calculate the angles at which shock waves propagate 


from supersonic aircraft. The generation of this waves can be seen in Fig. 2.9 and 2.11 


for a few different values of k= v/u. 





Figure A.| 


So we can say that this tail only happens when u < v. It starts to develop when 
u=v, like a shock wave, which is What we see in Fie) 2913; 


Let r be the length of the tail. 1t 1s easily computed: 
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Figure A.2 


We have that 
sin(0’/2) = u/v = R/r, so r= R(v/u). 


fine (t ) to dévelop this tailis t =r / v. 
The area of the tail will be, 2Area(ASO) - Area(AS’B) (Fig A.2). And 


Area(ASQ) = (1/2)R 1 
= (1/2)RV(? +R?) 
= (1/2)RSR?+ (Riva)? 
= (1/2)R*./1+(v/u)* 


Area(AS‘B) = a@R? and a@ = arccos(R/r). After substitution we can find the 


area of the target free tail. 


ea = Rey lq) - Rarccos(u/¥v) 
= Raley l +(viu)* - arccos(u/v)}. 


oF 


APPENDIAN B 
FUNCTIONS 


All functions tin this Appendix are used to calculate the shielded angle at any 
point in the area. In order to run these programs It 1s enough to enter target speed (u), 
searcher speed (v) and searcher radius (R) as a vector in function DENS Tie 
function runs all the other necessary functions as a subfunctions. The output density 


matrix can be found in variable RH. 
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Function ANGLE uses target speed(u), searcher speed(v) and _ searcher 
radius (R) as input and generates TIVV (8°) and T2VV (0°5). 

TRANSI TRANS2 and RRHO are used when the target is faster thanmeme 
searcher. Function TRANSI uses 0°,(TIVV) as input and transforms it to 8, (T1). 

Function TRANS2 uses 0°, (T2VV) as input and transforms it to 0, (Teas 
Function RRHO uses 9, (TI) and 0, (T2) as input to calculate densities p (RH) and 
puts Zeroes for densities on the detection cia. 

TRANI TRAN2 and RHO are used when the target is slower than the Seaneim 
Function TRAN uses 8° (TIVV) as input and transforms it to 0, (11). 

Function TRAN? uses 0’, (T2VY) as input and transforms it to oh). 
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Figure B.4 
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Function RHO uses 0, (Tl) and 0, (T2) as input to calculate densities p (RH) and 


puts zeroes for densities on the detection disk. 
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Figure B.6 





— 
—N 
By 
oc 
— o™~ 
N co~ 
By N 
oa x 
W S pe 
= ! om 
—N o™~ a 
BY re = 
oa o™ vi 
A + © 
: S 8 
Ds — x x —N ~ 
NTN on mine 
> ah ors Oo oO 
Sa N a ie = es van 
7 x Mm at N 
= A,A, = tha) 
eae ofe fe + © 
oe CONN ay | 
x NN wo 
pa ALA, ae ae 
PC + | N NN 
= dd Ay B® 
a AA,  & 
aati al n n 
a er coaro4es 
ay x x OOO 
see PaaS AOc 
4 
A 
A 
x 
r~ 
—N 
x 
=> 
A 
x 
a 
Je 


UxVx( (10K )x*2 ) 


2 


LONG bt tb ARYL sect 
*xXQdtdd IK AQUI 
MOD aS 


Xx k¥e~w OOOO NND 
m4 1X ®X NO, x 





—_— Y 

a a eee, ©, > 

oo Se a NO YRNRAWANNG 
eS Ch Ayey tt DO EEE 
mii iit I Ie 

Mmm trot NO AN MOT NOM OMOAN 


ANMANON OO AAA AA AAA AANNN 
[EL | SS | WES | | [SESE | WR | WN | | AT |W | W/E | | WW || | 


TETA2 PRIME ANGLE MATRIX WHICH IS 'T2VV' IN DENST FUNC. 


V DTRAN2 K:U:V;R: Pero eae St 
AS DEFINED IN DENST FUNCTION 


By 
v 
aa) 
N 


Figure B.7 
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APPENDIX C 
SINEULATION 


This simulation program is used to observe how the total area searched changes 


fea (Uuncion Ol U. Vv and t. 


PyoeeRCH JOB (4855,9999),'ROCKOWER ' ,CLASS=J 

V7 BAEC FORTVCG 

Py SORT ISYSIN DD * 

C SUBMARINE SEARCH AND DETECTION SIMULATION 
© PROE. &.. “—KROCKOWER WONTER 1967 
@ 


DIMENSION V0(7) OE NDET(30), XAXIS(0:30) 
REAL NSURY (0:30), DET (10 


Dea VOY ec a5 owe 2 a 2S 0 Bas 58) le Seay & 
U=1 

R= 1 

Te SOG 

YMAX = TMAX*U +R 

YMIN = -YMAX 


¢ 


UO(1) = .50331 
C DO TEN REPLICATIONS OF THE SIMULATION 

Dorda4 Re = 1,10 

ISEED = 337711*(U0(1) +.5) 

CALL LRND(ISEED, UO, 15000, 1, 0) 

TWPI = 2.%3.14159 

PemiOOeie= 1) 7 


= 770(T) 
K = 0 
BOme sek = 17.30 
NDET(K) = 0.0 
3)3 CONTINUE 


XMAK = (U+V )*TMAX + R 

XMIN = -R - MAX (0., (U-V ) )*TMAX 
AREA = (YMAX- 2urN) * (XMAX-XMIN) 
NUMITS= IFIX(AREA/3.) 
NSURV(0)=NUMITS 


DO 200 NUMIGT = 1, NUMITS 
M3 2 0. 

5 = al 
XT = XMIN + ae - IN} 200 (See 
YT = YMIN + *UO( 3*K-=-1) 


YMAX - YMIN 
Distoy= XTA*2 + YTA*2 

IF ( DISTO .LE. R**2) GOTO § 
THETA = TWPI*UO(3*K-2) 


UX = U*COS eae 
UY = U*SIN(THETA 
Die eR 104) (a VY) 


POes00.l- — 0 TMAK, DI 

XS = XS + V*DT 

XT = KT + UX*DT 

YT = YT + UY*DT 

Das i (XT-XS)**2 + YT**2 

1E( DVSieebe. R722) GOTO 10 

IF( (XT .GE. XMAX) . ORwetC. teehee XMIN) 2.ORe. (YT -.GE. YMAX ) 

/ OR. (YT. LIE, « YMIN) sOR, (DIST eGilk, DISTO ) GOTO 200 

PISTO = DIST 

300 CONTINUE 
GOTO 200 

10 IT = IFIX( T + .9999) 


(42 
ioe) 


200 


500 


100 
444 


/* 
// 


NDET( 1. )==— PETC eee 
CONTINUE 
DO 500 J = 1 


RAS 
CONTINUE 
eaCES (10) ate 
Dane = NUMITS-NSURV (30) 
WRITE(6,*) I, DET(I) 
CONTINUE 

WRITE(6,*) ISEED 
CONTINUE 

STOP 

END 


ee, 
NSURV(J-1) - NDET(J) 
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